Magnetization reversal in spin patterns with complex geometry 
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We study field-driven dynamics of spins with antiferromagnetic interaction along the links of a 
complex substrate geometry, which is modeled by graphs of a controlled connectivity distribution. 
The magnetization reversal occurs in avalanches of spin flips, which are pinned by the topological 
constraints of the underlying graph. The hysteresis loop and avalanche sizes are analyzed and 
classified in terms of graph's connectivity and clustering. The results are relevant for magnets with 
a hierarchical spatial inhomogeneity and for design of nanoscale magnetic devices. 
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Introduction. Reversal processes are of great impor- 
tance for technological applications, e.g., in information 
processing and memory devices. Often hysteresis curves 
with particular properties are required. In this respect 
novel "artificial solids" , arrays of nanoscale magnetic par- 
ticles 0, quantum cellular automata and integrated 
functional nanosystems Q offer challenging possibilities, 
yet to be understood. Currently methods are being de- 
veloped for patterning, defining and measurements of the 
magnetic properties at nanometer scale Q. Much less 
attention heve been devoted to the theoretical study of 
these systems 

In classical disordered ferromagnets 0, ferroelectrics 
and systems with structural transformations || the 
reversal processes are accompanied by avalanches, which 
are directly related with the motion and pinning of 
the domain walls. The distribution of avalanches and 
the time-series of the magnetization bursts (generalized 
Barkhausen noise) exhibits universal scaling features, 
which can be used for the diagnostics of the underlying 
domain structure. The pinning of the domain walls is vi- 
sualized as an effect of random disorder on a perfect crys- 
tal lattice, and commonly modeled by quenched random 
fields . Despite of large theoretical efforts 9] , the exact 
role of disorder in the emergent universality of noise has 
not yet been fully understood. In general, the stronger 
pinning implies smaller avalanches and reduced coercive 
fields (slimmer hysteresis loops). In the exact solution 
|To| in the case of Bethe lattice with large coordination 
number and weak disorder, a finite jump in the magneti- 
zation persists when the system size N — > oo, suggesting 
true criticality (absence of a cut-off) in the avalanche dis- 
tribution integrated along the hysteresis loop. 

Variety of domain forms may be nucleated in ferro- 
electrics with low symmetry and long-range elastic forces 
[ill ]. On the other hand, the patterned nanoscale mag- 
netic structures can be assembled in various geometries 
unrestricted by a crystal symmetry Q. The magneto- 
static interactions between neighbour nanoparticles and 
their shape anisotropy allow both ferromagnetic and an- 



tiferromagnetic coupling within nanoarrays [l2]. The 
full impact of the assembly processes and the emer- 
gent geometries on physical properties of the integrated 
nanoscale devices calls for theoretical analysis. 

The stochastic processes on complex networks at- 
tracted much attention recently I n general, the 
properties of the process depend on the complexity of the 
underlying network. Compared to more familiar Bravais' 
lattices, the network geometry makes severe constraints 
to the dynamics by restricting the interaction pathways 
along the inhomogeneous structure of links. The network 
inhomogeneity proves as beneficial in the case of trans- 
port processes Il4l . The spin dynamics on networks has 
been less studied |l5Hl(| . 

In this Letter we study dynamics of spins attached to 
nodes of a complex network and slowly driven by ramp- 
ing of the uniform external field. We assume antiferro- 
magnetic interactions between spins along the links and 
observe avalanches of spin flips due to topological inho- 
mogeneity of the network. An avalanche represents a 
fraction of spins that are reversed at current field value 
in order to minimize the energy within locally available 
geometry. We find numerically evidence of the hysteresis 
loop criticality which can be explained in terms of the 
structural properties of the network. 

Structure. Our model of two-state spins attached to 
nodes of a complex network can be regarded as a model 
of an integrated nanosystem with a nontrivial architec- 
ture. The structure is grown by systematic addition and 
linking of nodes 0|. At each growth step i a node is 
added and linked to M nodes selected among i — 1 pre- 
existing nodes. Selection of a target node, k, is given by 
a specified probability p(k,i). For p(k,i) = — 1), i.e., 
independent of properties of target nodes, the emergent 
network is known to have an exponential degree distri- 
bution. In contrast, the preferential linking, where the 
selection probability p(k, i) depends on number of already 
acquired links p( k, i) , are shown to lead to a scale- free de- 
gree distribution [ljj, Hj, . Here we use the probability 
p(k, i) = [l+q(k, i)/M]/2i to grow scale-free networks for 
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the purpose of this work. We restrict our discussion to 
the simplest complex networks for the following reasons: 
(i) Their structure is well understood 0, 0, 0] , in par- 
ticular, there are no hidden topological properties and 
link correlations [H; (ii) Varying a single structural pa- 
rameter M we tune the graph clustering property. Both 
scale- free network (SFN) and exponential network (EXN) 
have a tree-like structure when M = 1 (an example is 
shown in Fig. For M > 1 cycles (disregarding the 
direction of links) can appear, which strongly depend on 
the linking rule. For instance, for M = 5, we find the 
average clustering coefficient as 0.2803 in SFN, and 
0.0553 in EXN for network size of N = 10 3 nodes. 

The network sparseness and structural inhomogencity 
are the feature which can affect the spin dynamics. In the 
two network types introduced above, the profiles of the 
local connectivity (number of links attached to a node) 
are shown in the inset to Fig. [21 where i denotes the order 
of addition of the node to the network. In the case of SFN 
of size N the profile is given by the power-law dependence 
{q(i, N)) oc (N/i) 7 , with 7 = 1/2 related to the emer gent 
power-law degree distribution P(q) oc q~ 3 ( see [l5j, [l8j 
for more general scaling relations). In the case of EXN we 
have (q(i,N)) oc \og(N/i). When M > 1 the number of 
the elementary triangles or clustering (A(i)) at a node i 
(i.e., number of direct links among nearst-ncighbors of i) 
also varies through the network. The clustering profiles 
for M = 5 in SFN and EXN, exhibiting power-law tails, 
are also shown in the inset to Fig. [2] 

Spin Dynamics. A spin, S = ±1, is associated with 
each node on the network. The spins interact along the 
links connecting neighbouring nodes. We assume uni- 
form anti-ferromagnetic interaction J = — 1 and uniform 
external field H according to the Hamiltonian 

H = -J2 = -J2 JCijSiSj -H^Si- (1) 

The sum over pairs of nodes is restricted to the locally 
available links, which are given by the positive elements 
Cij of the adjacency matrix C of the graph. The inter- 
action along these links is symmetrical. Therefore, the 
local field hi, which is defined in Eq. varies with the 
local connectivity profile {q(i, N)) of the graph. 

We start from a negative saturation state in the field 
H = --ff max , where iJ max = QWx + 6 is determined by 
the largest connectivity <? ma x (at hub node) for each net- 
work realization, and S is a small shift which prevents 
hi = 0. Slowly increasing external field is outbalanced at 
some node where the interaction with other spins is the 
strongest. Flipping the spin at that node may launch a 
cascade of flips along the connections of the graph. The 
avalanche of flips is stopped when the energy is mini- 
mized at current field value. Then the field is increased 
again. This produces the step-wise increase of the mag- 
netization with time and with field, as shown in Fig. [3J 
in a manner known as Barkhauscn effect in systems with 




FIG. 1: (Color online) Scale-free tree with spins attached to 
nodes. Different collors correspond to spin-up (red) and spin- 
down (yellow) orientation when the field reaches H — + S. 
The snap-shot is taken during an avalanche. 

domain structure. Compared to the classical Barkhausen 
avalanches in driven disordered systems 0, H 01 ; a P ar t 
from the absence of disorder, the following important dif- 
ferences will be pointed out: (i) Only integer values of 
the local fields hi occur, making the lowest driving rate 
AH = 1 finite (ii) With zero temperature dynamics 
that we use, back flips are possible due to the antifer- 
romagnctic coupling; (iii) Disconnected avalanches may 
appear. We consider number of flips before a cascade 
stops as a measure of the avalanche size s. The net mag- 
netization changes may differ from the number of flips 
during one driving event. We define a time step as one 
update of the whole network. In the presence of the an- 
tiferromagnetic interactions the sequential updates are 
suitable for all M values. The distributions of avalanche 
sizes recorded along the ascending branch of hysteresis 
for the two network topologies and different clustering 
parameter M are shown in the main Fig. [21 

Hysteresis- Loop Criticality. We first discuss in detail 
the case M = 1, corresponding to a tree-like structures. 
As displayed in Fig. [21 the distributions of avalanche sizes 
exhibit a typical oscillatory behavior superimposed on an 
underlying power-law dependences, with different slopes 
in SFN and EXN. The characteristic peaks correspond to 
collective spin flips on a given depth layer of the tree 
(cf. FigHJ. In both cases the depth of a graph is pro- 
portional to D oc log(iV) (with different proportionality 
constants), however, the layers are more densely popu- 
lated in the EXN than in the SFN. For instance, the 
number of spins which are directly linked to the spin at 
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A V = SFN 




FIG. 2: Probability distribution of avalanche sizes s for differ- 
ent clustering parameter M =1, 5, 20. Two sets of curves cor- 
respond to scale-free and exponential networks with N = 10 3 
nodes, averaged over N s = 10 3 network realizations. Inset: 
Profile of the local connectivity (q(i)) (full lines) and num- 
ber of triangles A(i) (symbols) for the two network types, 
N = 3000, N s = 10. Fits: dashed lines (q(i)) oc log(10 4 /i), 
(A(i)> = 70(1 + i/10) -112 ; solid lines slopes: -0.5 and -1. 

the hub node (first layer) in SFN is sixfold larger than 
in the EXN of the same size (cf. inset to Fig. |2J . The 
population differences along layers roughly correspond to 
the relative positions of peaks in Fig. [5J 

Typically, spin at the hub is reversed early, then the 
structure around the hub remains stable until the field 
reaches very high values. This is compatible with the 
time picture, (Fig. Efb) shows M = 5 case), where the 
magnetization increase in the SFN takes much longer 
time, compared with the EXN. A snap shot of the spin-up 
and spin-down populations corresponding to the middle 
of the hysteresis loop in SFN is displayed in Fig. ^ Broad 
network connectivity distribution does not allow defini- 
tion of magnetic sublattices, such as familiar on Bravais' 
crystal structures. Moreover, the evolution of the den- 
sity of domain walls has a characteristic profile for each 
network structure, as shown in Fig. Ofb). This results 
in the hysteresis loop of particular properties. Fig. EJa) 
show the hysteresis loop for the SFN, which is dominated 
by the presence of the hub: Both, the extended narrow 
tails, and the discontinuity jump at the critical field, are 
attributed to the large connectivity of the hub node. 

With increased clustering of networks, i.e., when the 
parameter M > 1, closed loops of interacting spins ap- 
pear. The nearest- neighbor triangles have largest im- 



pact to the spin dynamics, being incompatible with the 
antifcrromagnctic interactions. This leads to frustration 
of spins, a familiar concept in the theory of spin-glasses 
|2l[. However, the observed inhomogeneous distribution 
of triangles on the scale-free graph, as shown in the inset 
to Fig. [21 affects the hysteresis loop criticality in differ- 
ent manner, compared with the case of spin-glasses |22j . 
In particular, the critical field H c of the magnetization 
reversal increases with the clustering on the scale-free 
network and a finite magnetization discontinuity persists 
even for very large clustering (see Fig. Efa)). On the 
other hand, the absence of the hub nodes in a mild loga- 
rithmic profile of the EXN, results in a narrow hysteresis 
loop and the absence of the large discontinuity for all 
M > 1, more similar to the case of Bravais' solids |'22l. 

Increased clustering results in stronger pinning |l6| . 
Consequently, the cut-off decreases in the avalanche size 
distribution, with gradual smearing of the oscillations. 
The distributions approach the form that can be fitted 
by a stretched exponential function for M = 20. The 
qualitatively same behavior with different parameters is 
found in both types of network structures (cf. Fig. |5J). 
The distributions of avalanche durations are broad with 
increased cut-offs for large clustering. However, owing 
to small depth ( "small- world" feature) of the networks, 
the overall range of the distributions is short (about one 
decade) to determine a precise functional form |2(| ■ 

In order to understand the nature of the hysteresis 
loop criticality in driven disordered ferromagnets, re- 
cently several exact results on Bethe lattice 01 an d other 
sophisticated theoretical approaches were applied [2j| . It 
is likely that these approaches, based on the analogy of 
the Barkhausen avalanches with the stochastic branch- 
ing processes and bootstrap percolation, can be useful 
for the study of the hysteresis curves on networks, with 
a necessary adjustments that will take into account two 
essentially different points: (i) antiferromagnetic interac- 
tions and (ii) power-law connectivity profile. In particu- 
lar, based on our numerical results in Fig. we expect 
that true criticality can be proved in the case of tree 
structures (M = 1) with the scaling exponents of the 
avalanche sizes expressible in terms of the connectivity 
profile 7 as r s = 1 +7 = 1.5 for the scale-free, and r s = 1 
for the logarithmic profile in the exponential tree. 

Conclusions. We have shown that topological con- 
straints due to complex geometry of network and anti- 
ferromagnetically coupled spins leads to the avalanche- 
like magnetization reversal and hystcrcsis-loop criticality 
without disorder. The broad distribution of the network 
connectivity and the presence of the hubs in the scale-free 
structures leads to unexpected features of the magneti- 
zation reversal: increased clustering promotes stronger 
pinning and shorter avalanches, but at the same time 
increases the critical field and preserves a finite disconti- 
nuity in the magnetization. These features may require 
renewed theoretical concepts. In practical terms, our re- 
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FIG. 3: (a) Hysteresis loop for the scale-free network of different clustering parameter M = 1 (tree), 5, and 20. (b) Time 
dependence of the magnetization, number of flips (Barkhausen noise), and domain- wall density (broken lines, shifted by —1 for 
better view), for two types of network structures and for M = 5. 



suits suggest that the control of the pattern geometry im- 
plies the control of the hysteresis curve properties within 
an enlarged range of parameters. Our simplified model 
with "free graphs" may initiate study of more realistic 
models and practical analysis of integrated nanosystems 
with rich structural and magnetic characteristics. 
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